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ABSTRACT 


The generation of spontaneous magnetic fields in a 
plasma produced by a 5.0 Joule, 22 nanosecond laser pulse 
and the generation of the reversed magnetic fields after 
several hundred nanoseconds due to the interaction between 
the laser produced plasma and a background plasma was 
investigated. 

The generation of the reversed magnetic fields coincides 
with the steepening of the expanding plasma front. This 
steepening is caused by a collisional shock wave which is 
formed during the interaction of the expanding laser- 
produced plasma with the photoionized background gas. The 
analysis of the source term for magnetic field production 
indicates that the generation of the reverse fields is due 
to an axial temperature gradient and a radial density 
gradient. This axial temperature gradient 1s a result of 
the heating of the plasma in the collisional shock wave. 

The width of the shock was measured to be 0.2 cm. The 
magnitude of the maximum reverse field was measured at 5.0 


worr (H.) background pressure to be 150 Gauss. 
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toe ERODUCTION 


A dense, energetic plasma can be produced by focusing 
the pulse from a powerful laser onto a solid target. Asso- 
ciated with this laser-produced plasma (hereafter referred 
to aS laser plasma) spontaneously generated magnetic fields 
liave been reported and investigated [Refs. 22,24]. These 
investigations imply that spontaneously generated magnetic 
fields arise generally during the production of a laser 
plasma. The strength of these fields may vary with target 
geometry [Ref. 31], background gas pressure [Ref. 3] and 
amount of energy delivered to the target by the laser. [In 
a computer Simulation experiment, Widner [Ref. 30] has 
observed a spontaneous field intensity of the order of l 
Megagauss. 

Stamper et. al. [Ref. 25] have reported spontaneous 
magnetic fields, of the order of 0.6 Megagauss, generated 
in the plasma near the focus of a high-irradiance Nd-laser 
beam. They used a Raman-shifted, three channel Faraday- 
rotation diagnostic system and observed that the spontan- 
eously produced magnetic fields depend on a variety of 
experimental parameters, such as target material and 
geometry, laser irradiance and polarization, focal position, 
and prepulse effect. 

The influence of magnetic fields on laser plasmas is an 


important topic in laser-fusion research. The presence of 


i 





such fields in a plasma has effects on electronic heat 
conduction and the generation of plasma instabilities. 

After laser shut-off, generation of reverse fields has 
been observed at the front of the expanding plasma for 
ambient pressures above a certain critical value [Refs. 3 
ema 21). 

This project contains the study of field reversal in 
an expanding laser plasma. 

The spontaneous and reversed fields were detected using 
magnetic probes (inductive loops). The probe coils con- 
Sisted of five turns of #40 FormVar copper wire. 

The plasma density profiles were diagnosed with elec- 
tric double probes. The double probes were maintained at 
a fixed bias of sufficient magnitude (20 volts) to insure 
collection of saturation lon current. Since the expansion 
velocity of the laser plasma was greater than the ion 
thermal velocity, the probe current was proportional to 
the product of the ion density and the plasma expansion 
velocity. The expansion velocity was obtainable from the 
probe signals so that the relative density of the plasma 
mould be calculated. 

This thesis is divided into 7 more sections. Section II 
contains the previous experimental work connected with the 
Field reversal. Section III discusses the basic theory of 
spontaneous magnetic fields in laser produced plasmas, field 


eeéversal and shock waves ina plasma. Section IV contains 


Ibs 





the details of the experimental methods. The experimental 
results are presented in Section V. In Section VI, the 
experimental results are compared with the theory. Section 
VII contains the summary and conclusions. Recommendations 


for further study are contained in Section VIII. 
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ite ee OOS RE ORARCH 


Peeing  EMestatter pase Of 1970, aS part of a study of 
the interaction between a laser plasma and the ambient 
background plasma, Stamper et. al. [Ref. 24] discovered 
that magnetic fields were spontaneously generated in the 
laser plasma. They explained the generation of these 
fields in terms of thermoelectric currents associated with 
large temperature and density gradients near the target. 

The spontaneous magnetic fields of the order of a kilo- 
gauss, were observed with magnetic probes. After this dis- 
covery, McKee [Ref. 18] made extenSive studies of such 
fields at the Naval Postgraduate School. During McKee's 
investigation the following observations were made: 

(1) The maximum fields detected at a fixed position 
showed a systematic dependence on the ambient gas pressure. 

(2) The fields were symmetric about the target normal. 

(3) After several hundred nanoseconds, a relatively 
small magnitude field in a direction opposite to the direc- 
tion of the initial field appeared behind the initial 
fields. This reverse field appeared sooner in time for a 
Mmmeeground pressure of 250 m Torr than for 0.1 or 5 m Torr. 

McKee developed a model to explain spontaneous magnetic 
fields. According to McKee's model: If a background 
plasma is present, the expanding laser plasma will couple 


to it. If the density of the background plasma is high 


IS 





enough so that the mean free path for binary collisions 
between the laser plasma ions and the background plasma 
lons is smaller than the dimension of the laser plasma, 
the coupling will be collisional in nature. At lower 
background pressures collisionless coupling will arise 
because the presence of the spontaneously generated magnetic 
field at the front of the laser plasma will lead to coupling 
Via the ion-ion two stream inStability in the presence of 
a magnetic field. As a result of the coupling ambient 
plasma will be swept-up at the front of the expanding 
plasma and a density shell will form. The shell contains 
geeeensity gradient in the z-direction (see Figure 1); then 
field generation can occur in the density shell if a radial 
temperature gradient arises (see Section III for details). 

After McKee, Bird developed another model (interaction 
model) to explain the pressure dependence of spontaneous 
magnetic fields. This model attributes the increase in the 
magnitude of the initial spontaneous magnetic fields with 
increasing background gas pressure, to the interaction which 
occurs between the laser plasma and the ambient background 
plasma. The interaction increases the axial density gradi- 
ents of the front of the expanding plasma and this, in turn, 
increases the rate at which spontaneous magnetic fields 
are generated in the interaction region. 

During Bird's investigation no density shell was detected. 


Also the shell model could not explain the field reversal. 
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The following observations on reverse field generation 
have been reported by Bird: 

(1) Field reversal is coincident with a strong inter- 
action between the laser plasma and the background plasma 
aS is evidenced by a steepening and slowing of the front 
of the expanding plasma. 

(2) The onset of field reversal occurs at the front 
of the expanding laser plasma and it occurs earlier in 
Space and time at higher pressures. 

(3) Reverse field can be produced long after laser 
shut-off if the laser plasma encounters a barrier. 

(4) Reverse field iS initially produced at the front 
of the laser plasma streaming into a backgroundof H, at 
pressures above 450 m Torr before any Significant pile-up 
of the laser plasma has occurred. This narrow Pon monente 
of reverse field appeared to be generated in the swept-up 
hydrogen background plasma and it was followed by the 
appearance of a broad component of reverse field as the 
plasma front steepens and slows down due to the interaction 
between the laser plasma piston and the background. 

These results suggest that the generation of the reverse 
Magnetic fields can be explained by formation of a shock 
wave driven by the expanding laser plasma piston. Such a 
Shock wave compresses and heats the plasma. The formation 
and the effects of the shock waves in a plasma are discussed 


foecne LOllowing section. 
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III. THEORY 


A. LASER PRODUCED PLASMAS 

The production of plasma by irradiation of a solid 
target was first treated theoretically by Basov and Krokhin 
feet. 1]. 

The plasma is created by heating, vaporization and 
1onization of the target material. The overdense plasma 
cs > no)” at the surface of the target expands until the 
Pmeceron density in the front drops below the critical 
density Nos Strong absorption occurs when the density 1s 
near the critical density ae and the plasma is heated to a 
high temperature and the pressure gradient accelerates the 
plasma. After the laser shut-off, the plasma front expands 
at a constant velocity. Thus, three distinct states can be 
observed during the course of the expansion of the laser 
plasma. 


When the expansion is into an ambient background gas 


rather than into vacuum, the last stage looks quite different. 


tthe critical density n_ for the laser plasma is the 
density at which the electr6én See 2 1s equal to 
the frequency of laser radiation M& /e*). At elec- 
tron densities above this . value, ths laser plasma 
reflects most of the incident radiation. pane emirerca | 
@emsity for 1.06 um radiation is about 10 em 
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At this stage, for a background gas pressure greater than 
0.2 Torr, Bobin et. al. [Ref. 2] observed a luminous sheet 
ahead of the expanding plasma. This luminous sheet lasted 
more than 150 nsec after the laser pulse peak value. Its 
velocity decreased with time. The authors attributed the 
Origin of the luminous sheet to a shock wave driven by the 
plasma acting aS a piston, because their observed time evolu- 
tion of the luminous front agreed both qualitatively and 
Quantitatively with a simple shock wave model. 

A one-dimensional, simplified model treating the expan- 


sion of the laser plasma into vacuum has been given by 


Ready (Ref. 20]. In this treatment it is assumed that: 
(1) The plasma expands in the z-direction. 
(2 Incident radiation 1S absorbed via the inverse 


Bremsstrahlung process (free-free absorption by electrons 
in the laser plasma. 


eesnewmescget 1s inicially a codd fullyeaionized 


plasma. 
(4) Radiation losses can be neglected. 
(5) The electrons and ions have a common temperature. 
(6) Temperature is independent of position. 
2 


The equipartition time for transfer of absorbed energy 
from electrons to ions can be approximated by 


eg _ aL By ae —3 -l 
Peo) Te (°K) [In fem ~)2] 
where m. is the ion maSs and m_ is the Proton mass. For 
n = 10*+ ana Te = 100 ev, a ST aascc - 


aS, 





(7) The plasma expansion velocity is of the form 


pySically, the pressure gradients within the laser plasma 
are driving the expansion. The form of the pressure 


Variation can be approximated by 


where oe is the pressure near the surface of the target, 
and Z is the distance of the plasma front from the target 
surface. 

The temperature of the plasma can be determined by 

‘ 

equating the rate at which the laser beam deposits energy 
in the plasma, to the rate of change of plasma thermal 
energy plus the rate at which work W is done by the plasma. 
For one dinensional plasma of area A, expanding in the z- 


Seeection this gives: 





ic Wee dP az 
Pca eode let <? S55 Ge oo 
From the plasma equation of state: 
Pai Ge nS) ee Oy, 


mmemcan Obtain the variations of electron and ion densities: 
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where ae and Nog are the lon and electron densities at the 
target surface. The work done by the plasma in the expan- 


sion increases the plasma kinetic energy so that: 


dP CZ ei, Cl Cs 
32° ae OO) Ge la (nymytngm,) (Fe) I so) 


Beeantegrating Equations (1) and (4) from zero to Z and 


the result, 
“a 
No SR if eth @ Sicke  Npa/ 
S Z, S : 
0 


where N is the total number of ions or electrons in the 
plasma and Ns 1s the corresponding number density at the 


target surface, Ready obtained the equations 


(N.m.+N m )Z Pa 
11 eee 





dz 
ee eS Se ee (5) 
6K(N, *N_) ae 
a 
q a 5 DA ON a (6) 
ale Le eae e @e 


These two coupled equations were solved numerically by 
McKee for the conditions of this experiment. The maximum 
temperature obtained was of the order of 100 eV and this 
maximum occurred a few nanoseconds after the arrival of the 


maximum laser intensity at the target. 
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In the above treatment, 1t was assumed that the 
expansion of the plasma is into a vacuum. The presence 
of a background gas (pressures above 1m Torr) influences 


the expansion dynamics. 


iB . INTERACTION BETWEEN COUNTERSTREAMING PLASMAS 

The optical radiation from a dense, hot laser plasma 
1s known to be rich in the far ultraviolet and near x-ray 
regions of spectrum. Line radiation and Bremsstrahlung 
radiation from a plasma of T = 100 eV is energetically 
capable of ionizing an ambient background gas. 

Koopman [Ref. 16] has investigated the photoionization 
produced from a plane copper target by an 8-Joule, 20 nsec 
laser pulse at 6943 A° from a Ruby laser. He found that 
the ambient gas was almost completely ionized to distances 
of the order of a centimeter for a background gas pressure 
of the order of 0.2 Torr and at large radii the degree of 
photoionization fell as a Total optical energy emitted 
between 1000-2000 A° was eStimated as approximately 170.5 
Joule. Therefore, when a laser plasma expands into a 
photoionized ambient background gas, plasma-plasma inter- 
action must be considered on the expansion dynamics. 

According to Dean [Ref. 7] during the period of the 
laser irradiation of the target (early time), the dynamics 
of the laser plasma expansion corresponds to the propagation 
of a radiation-driven detonation wave. After the end of 


the laser pulse (late time) the plasma expands as a blast 
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wave. This model assumes that at early times of the 
expansion phase, there is momentum transfer between the 
expanding front and the ambient background plasma and that 
as the front propagates into the background plasma, further 
energy from the laser is absorbed behind the front in which 
the density is near the critical density and is transported 
to the front by electron heat conduction. After the end 

of the laser pulse, all the momentum in the driving laser 
plasma piston is assumed transferred to the swept-up back- 
ground plasma which propagates as a blast wave. 

Hall [Ref. 12] investigated the dynamics of the expan- 
Sion of the laser plasma into an ambient background of 
eagon. During his investigation, a 0.5 Joule, 20 nsec Ruby 
laser was used. The target was a Tantalum slab. Using an 
image-converter camera and double probes for diagnostics, 
he found that the late time expansion dynamics of the laser 
plasma could be explained by blast wave theory. 

Koopman [Ref. 16] has suggested that the observed momen- 
BemecGOupling is the result of multi-coulomb collisions 
occurring between particles of the two plasmas. If a strong 
momentum transfer interaction occurs, one would expect 
"Snowplowing"” or displacement of the ambient plasma by the 
leading edge of the laser plasma piston, with a corres- 
ponding pile-up of piston plasma behind the interface between 
the laser plasma and the background plasma. The ambient 


and piston plasmas would be intermixed in the region of 


is & 





interface, that is, the model is one in which the background 
plasma is "captured" at the front of the expanding plasma, 
increasing the total mass of the streaming plasma and, by 
conservation of momentum decreasing the expansion velocity. 
A one dimensional expansion model can be justified on 
the basis of the following considerations. At the target 
surface the laser beam diameter is about 1.5 mm and the 
laser radiation is absorbed within a distance of the order 
@eeo.01-O0.1 mm from the target [Ref. 20]. Therefore, the 
ablated target material initially will have larger dimen- 
ees In the radial direction than in the axial direction. 
Since the expansion of the laser plasma is driven by pressure 


gradients; 


Pp Pp 
Cre — and [ules = 
Z ie 
from the above discussion, a >> ae therefore 
mere). {| >> ee whiccenon and 6, are characteristic 
radial and axial scale distances. From the above inequality, 


a one dimensional laser plasma expansion can be approximated. 
When the laser plasma expands into an ambient background 
plasma, the axial position of the expanding front will be 
governed by the interaction between the two plasmas. An 
empirical expression for the late time (after laser shut- 
off) dependence of the axial position of the plasma front 


On time and on the ambient (nitrogen) plasma density, pee 
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moegivyen by Bird {Ref. 3]: 


+. 


Po 


0.10.5 


48) a=) 


Gee LNITIAL PRODUCTION OF SPONTANEOUS MAGNETIC FIELDS 
In the presence of electric and magnetic fields, the 
equation of motion for the electron fluid is given by 


generalized Ohm's law. 


nm [——= + (V_-V)V_] = oe ae > ee (3) 
where Vi (r,t) fomeiemrruicevyolOcicy Of electrons, E(r,t) 
and B(r,t) are the averaged electric and magnetic fields, 
i is the kinetic stress tensor and ee represents the 
collisional momentum transfer between electrons and ions 
per unit volume and unit time. 

If one assumes that the electron-electron collision time 
Pe. is much less than the characteristic time for field 
emeeuccion which is of the order of duration of the laser 
paise, then ier reduces to the electron pressure gradient, 
VP. Pemtie mMacrOoscople Llow velocity of the electron 
Fluid element is small compared to the random thermal 
velocity of the electrons then the ee) term can be 
neglected. For the aluminum laser plasma, the expansion 
Memecity OL the front is of the order of 10! cm/sec 


Peet. 5]. If one assumes an electron temperature of 100 eV, 


ZS 





then 


Vv. | 


Se Se OTL 
1/2 
[kT /m,] 


Therefore the above approximation is justified. 

If the frequencies of interest are much less than the 
electron plasma frequency, 1.e., w << Woe the electrons 
can be considered massless and react instantaneously to 
applied fields. 


Pa can be approximated by the resistive drag term, 
ee) 


mewiiecre GC is the conductivity of plasma. Finally, 





eq. (8) can be written 


Byeecaxing the curl of E and using the Maxwell equations, 








eq. (9) becomes 
-VxE = SB =~ Syxo+ 9 %(V_xB) + Vx UP (10) 
ce O e ; ne e 
Also, from the Maxwell equations, 
- l Oo: 
VxB = Hod ac a2 TE (LL) 


The ratio of the two terms on the right hand side of 


eq- (ll) is of the order 
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where as 1s the duration of the laser pulse ae 108” sec) 
Bereeal] situations of interest o > 1 mho. Therefore, 


eq. (12) becomes: 








Ee Ee ai kal 
oe ee lg tena El < La os ae i 
2 OT 4g T 
eA oe O O O 
oe) 
The above inequality shows that 
1 28 
32 dt 
<< L 
Hod 
i.e., the displacement current can be neglected in compari- 
Pemero the conduction current. In doing so, equation (11) 
becomes 
aS Wo J ies; 
meng eq. (13), the first term in eq. (10) becomes 
ah eee! , ® Se eee oe 
2 VxJ = Tach = i aluy B VB) i} 
O O 
(ks) 
Since 
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eee OM - Saxo = —. o* B (14) 
0 il Ie; 
O 
imae chird term in eq. (10) can be written: 
kT 
dk = iE _ k = 
ea’ Fo = V x ao Meet e! = V Sa a n = ile 
e e e 
KT. 
= V Sl eaerah st (dey, 
e 
as 
The density gradient is a vector; —— = A but the tempera- 
kT e 
ure 1S a scalar quantity; — =O echen, USinNg the vector 
identity; 
Vx(ZA) = @VxA + (VM)xA icleG)) 
Eq. (15) becomes: 
kT 
Vx(—~20Vn_) = Sur xLyn (17) 
noe @¢é 2a 


Meacting ali terms in their final form together, eq. (10) 


becomes 


ie ee es i 
= Tar B 4 Gy x2) - Bip) oe ae Gis) 


= 


Qo 


Qo 
a ee 


The first two terms on the right hand side of eq. (18) 
are the diffusion and convection terms. There is no external 
magnetic field applied, i.e., initially B = 0. Therefore 


the generation of an initial magnetic field requires that 
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the last term, the thermal source term (S) be non-zero 
due to non-parallel temperature and density gradients. 
The initial configuration of magnetic field will decay 


away in a diffusion time t [Ref. 14], 
(19) 


where Lis a length, characteristic of the spatial varia- 
tion of B. For times short compared to the diffusion time 
Or, in other words, when the conductivity 1s so large that 
the first term in eq. (18) can be neglecged, equation (18) 
can be reduced to: 


Yaak Vn 


Se Se 
Sees ea) x ) (20) 
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This equation expresses the fact that the rate at which the 
fields vary in a region, depends on the rate at which the 
fields are created in the region and the rate at which the 
md £iowS in or out of the region. 


For adiabatic expansion of the electrons, pv’ = const. 





M 
(ee eS 


omen kT 
e "e’m n 
ao 


const. (where M PS cheweocal electron 


i= 


mass in a given volume V) or KT const. and the 


source term vanishes en cm AV. =O) Lous, non— 
adiabatic conditions are necessary for the generation of 
magnetic fields by the thermal source term. 


If the laser energy is supplied to the electrons at a 


Sufficiently high rate, then the magnetic source term is 
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non-zero in the resulting non-adiabatic state. Stamper 
[Ref. 26] analyzed the problem by means of the electron 
heat equation. 

Gur dn 


el Se se dn 
con a SS aeke 


where W,, is the energy supplied to the electrons per unit 


‘a 
time. The first term on the right hand side of equation (21) 
ms. CT al ' where W. is the internal ener and the 
y-l dt dt ic SY 
second term on the right hand side of equation (21) 

dw. 
Kren = where W.. 18s the work done by the electrons. 


dt dt W 


meeording to Stamper, the ratio W,/W, is the measure of 
Bemeadlabatic conditions. Stamper approximated Wy by th 
Beeduct of the absorption coefficient yy Cone) and the 


Baddatilon intensity I (Watt/cm*) , while W,; was approximated 
NekT 


Ph 
one assumes k,, is due to inverse Bremsstrahlung absorption, 


by 





where oe is an average electron heating time. If 


in the underdense region, that gives: 
: =! _ 
ee “@il Ree en 
Tos NkT 


We el 








(22) 


where To 1s the electron-ion collision ie. i. cheeeelev= 


J: 


eee te? °K) for n. ~ 10°> em > 
el = 3 un 
emis 2) 
A 
and T. ero Mey. From footnote (2), 1 es se sec. a 
jme~ 8] x 10° 1° SCemeonemtnie Tatient:./T_.. = 9.25 x 10° ee 
el el eL. ~ m 
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index of refraction of plasma (n ~ 1). Actually near the 
cut-off, the index of refraction goes to zero, thus the 
ratio W/W >> 1 and the plasma appears non-adiabatic. 
During the time the plasma is absorbing the laser radiation, 
for the conditions of this investigation, the ratio W/W 
is of the order of one and the plasma appears non-adiabatic. 
em: if Wr, << We and Wr then the plasma appears adiabatic. 
The direction of magnetic field is determined by the 
geometry of the laser plasma. At early times, the laser 
plasma, as justified before, will expand in the axial direc- 
tion and create an axial density gradient. The temperature 
gradient is due to the finite radial extent of the laser 
beam and is in the negative radial direction. The cross 


product of these two gradients will generate azimuthal and 


axisymmetric magnetic fields [Fig. 6] as observed. 


Peer ie LD REVERSAL 

At later times, after cessation of the laser pulse, 
the gradients become smaller and field production becomes 
insignificant. However if the laser plasma expands into 
a background plasma of sufficient density, shock heating 
at the front will produce an axial temperature gradient 
[Ref. 21]. Since a radial density gradient always exists 
in the plume-shaped laser plasma, magnetic field in a 
direction opposite to the initial field will be generated. 


The equation describing this field is: 








) C28) 
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Schwirzke [{Ref. 21] and Bird [Ref. 3] observed this magnetic 
field reversal several hundred nanoseconds after laser shut- 
off. The production of reverse field at the front of the 


laser plasma requires that the ratio: 


ALA 


ere dle 


be greater than one. 


The gradients in this ratio can be approximated by: 


where the 6's are the characteristic gradient scale lengths. 


The ratio becomes Sean Ope ES 4 -. For a one dimensional 


AR 2T) iam 

adiabatic compresSion or shock compression in the z-direction, 
=) 0 << : - 

5 or oe and oon 5 em The large electron heat conduc 


tion provides that 5. > and thus the ratio becomes 


T 
greater than one lie., reverse field production is possible 
under these conditions. 

The production of reverse field appears to be a conse- 
quence of the interaction between the laser plasma and the 
ambient background plasma. The interaction results in 
slowing down and a steepening of the front of the laser 
Plasma (pile-up). 

Mauiation (23) predicts that the field reversal will occur 


sooner and grow to a larger magnitude as the background 


pressure increases. Bird [Ref. 5] justified this prediction 


a2 





and also observed that 

(1) The mechanism for production of reverse fields 
is "turned-on" sooner at higher background pressures. 

(2) Field reversal iS coincident with a strong inter- 
action between the laser plasma and the ambient background 
aS 1S evidenced by a steepening and slowing of the front 
of the expanding plasma. 

The electron temperature required to produce a given 
field reversal can easily be estimated by dimensional 


analysis of eq. (23). From that equation 





where Se and 5 are radial and axial characteristic scale 
lengths and t is the characteristic time for reverse 
magnetic filed production. Bird estimated the required 
temperature to be of the order of 1 eV for the parameters 
of his investigation which were similar to the ones of the 


present investigation. 


E. SoHOCK WAVES IN A PLASMA 

imeec large concentration of energy in a continuous 
medium is suddenly released, it will spread into the medium 
and at its forefront will be a shock wave. A shock wave is 
a Sudden but continuous change in the state of medium, heats, 


compresses and accelerates the medium through which it 
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paepagalLes. The shock wave may propagate through a 
plasma. 

The normal process of shock formation is the pressure 
dependence of wave velocity, so that a finite amplitude 
disturbance with increasing pressure steepens as the higher 
pressure portions move faster than the smaller ones in 
front. In order to obtain a steady profile wave, some 
mechanism must exist to limit the steepening. This limi- 
tation on steepening can be provided by either dissipation 
Or dispersion. A steady state shock occurs when a diffu- 
Sion process balances the steepening. Jump conditions are 
then obtainable from the conservation of mass, momentum 
and energy. In gas dynamics, the flows are properly des- 
cribed by the Navier-Stokes equations. Any shock wave 
present in the flow must therefore also be contained in the 
solutions of the Navier-Stokes equations, under appropriate 
maettal and boundary conditions [Ref. 3]. However if the 
thickness of the shock wave becomes so small as to be 
comparable to the mean free path, the Navier-Stokes equations 
become invalid in describing the flow in the shock. One 
must then use the Boltzmann equation or similar kinetic 
equations for the description of the shock structure. 

The basic features of the shock structure in a plasma 
are related to the slow character of the energy exchange 
between ions and electrons and to the high electron mobility, 


ema CFesult of which the electron heat conduction greatly 
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exceeds the ion heat conduction. Maxwellian distributions 
in the electron and ion fluids are established quite 
rapidly, in atime of the order of the time between parti- 
cle collisions [Ref. 33]. On the other hand the equili- 
bration of the temperatures of ion and electron fluids 
takes place much more slowly, because of the large differ- 
ence between the electron and ion masses. This relaxation 
process determines the shock front thickness in a plasma. 

A qualitative discussion of the consequences of the 
low rate of energy exchange between electrons and ions is 
fepsen by Zel'dovich and Railzer [Ref. 33]. In this simpli- 
Fied discussion it is assumed that: 

(1) electron heat conduction does not differ from ion 
heat conduction 

(2) ionization does not take place in the shock wave; 
the wave 1S propagated through a gas which is already 
Lonized 
In a reference frame moving with the wave a considerable 
part of the kinetic energy of the gas entering the com- 
pression shock, 1S irreversibly converted into heat through 
the action of ion viscous forces. The increase in lion 


temperature across the compression shock is of the order 


m. V 
aa 


k 


equal to the shock front velocity and m; 1s the mass of an 





AT. ~ where V is the velocity of the incident gas, 


ion. The thickness of the viscous shock is determined by 


the time between ion colliSions, Tai 1t is of the order of 
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an 1lon-ion collision mean free path Re ~Vienta where Wie 


is the thermal velocity of the ions in the compression 
Shock. During the compression time tT. the ions do not 
Succeed in transfering any appreciable thermal energy to 


the electrons, Since the characteristic exchange time 
E 


m. 
io. ~ pty 172. 
el mM, aL 


j=me T.. 
oe 


is very large compared to the ion collision 


The increase in electron temperature across the shock 
due to the conversion of the kinetic energy of the elec- 
trons into heat through the action of electron viscous 
memees 1S negligibly small. it is of the order of 


AT. = meus 
k 


The electron gas in the compression shock is heated by some 





, that is smaller by a factor of m/™; than AT, . 


other process. Since the dissipation of energy due to 
electron viscous force is very small, Zel'dovich and Raizer 
(Ref. 33] suggested that the compression and heating of 


the electron gas in the compression shock is adiabatic: 


Te2 etl 
= = a (25) 
el iL 
where we and T32 are the electron temperatures in State I 


and State II, Py and P, are the electron densities in State 


I and State II and y is the ratio of specific heats. For 


0 
@ strong shock wave, the ratio = Po wOrecene order of 4. Thus 
AL aL 
aoe Mees bor y = 1.4. Therefore in a strong shock wave 
el 


propagating through a plasma with equal ion and electron 
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temperatures, behind the shock there is considerable 
difference in the temperatures of the two gases. After a 
small portion has undergone the shock compression it begins 
to transfer thermal energy from the ions to the electrons 
and this leads to an equilibration of the temperatures in 

a time of the order of the energy exchange time Baa The 
final temperature Ty (= ee = T;) is determined by the general 


equations of conservation for shock fronts. The thickness 


of the relaxation region behind the shock in which plasma 


approaches equilibrium (T. = T; = Ty) is of the order of 
p 
Ax a aa where Y 5 oli —vV 


Jaffrin and Probstein [Ref. 15] investigated the structure 
of a one dimensional shock wave ina fully ionized plasma 
consisting of electrons and Singly ionized ions in the 
absence of any applied external electric or magnetic fields. 
The analysis of the shock structure is based on a two- 
mora Continuum approach in which they took as a "model" 
the Navier-Stokes equation, including the electric field and 
transfer terms. In so doing the overall width of the shock 
transition region turns out to be considerably larger than 
the characteristic mean free path. They found that 
(1) for mach number M < 1.12 (weak shock) the shock 
Drorile is monotonic, with the shock thickness of the order 


OL LaSe where Le, is the down stream (behind shock) mean free 


So 
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waed Of LOn-1oOm Collisions and € = Cn . The important 


dissipative mechanisms are electron ene eonaucelon and 
the energy exchange between ions and electrons. 

(2) M > 1.12 (moderate to strong shocks). There 1s 
Still a thermal shock layer whose approximate thickness 
is Loe. In this layer the electron temperature increases 
and approaches its downstream value while the lon tempera- 
ture and velocity changes are relatively small. However, 
imbedded in this first layer there exlsts a second inner 
shock whose thickness is of the order of the downstream 


1on-ion collision mean free path 2 In this second layer, 


5° 
the dissipations are due to the lon viscosity and thermal 
conduction. The heating and compression of ions occur 
mainly in this second shock layer. However, because of 
the ions' much heavier mass, the ions do not have time to 
exchange energy with the electrons in a few mean free paths 
and the electron temperature 1S essentially unaffected. 
The ion and electron temperatures then equilibrate in the 
downstream continuation of the thermal layer. The thick- 
ness of this downstream layer is of the order of Lao/e. 

(3) Debye length of the order of the mean free path 
(A, 2 ha). In this case spatial oscillations of wave- 
lengths which are of the order of the Debye length are 
present. The charge separation is no longer small and the 


electrical effects interact with the flow in the two 


— electric shock layers. 
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The structure of a collisional weak shock in 
fice presence Or magnetic Etteld which is in the plane of 
shock was studied by Hu [Ref. 13]. He found that 
the shock thickness depends on four parameters, the 
electron mean free path er the electron larmor radius rr. 
based on thermal speed, the electron penetration length 
yee. 
e 


r= on 17>” and the electron-ion mass ratio ao. Hu 
P (477 .e ) 


discovered two important features of the shock structure. 

(1) It is possible to have collisional shocks with 
oscillatory structure due to electron inertia. 

(2) It 1S possible to have collisionless shock with 
collisions whose thickness is smaller than the collision 
mean free path due to the complicated interactions between 
electromagnetic fields and the plasma. Often in interpreting 
experiments, shocks exhibiting these two features are con- 
Sidered as collision-free shocks. 

In general, the mean free path for collisions ina 
plasma determines the type of the shock. A shock layer, 
which depends on the interaction length, can be thinner 
than a mean free path. These shocks provide examples for 
Studying the mechanism of collisionless dissipation. A 
Similar definition of the collisionless property was given 
by Friedman [Ref. 10}. For a given shock structure, the 
dissipation provided by classical collisions (ordinary 
momentum transfer collisions) must be less than the dissi- 
Pee On which is required to form a steady shock. 

Mie £lmst conclusive evidence of the collisionless 


shocks come from astrophysics. The plasma flow from sun, 
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the solar wind, is highly supersonic when encountering the 
magnetic field of the earth, and a shock wave must be 
formed. However, the density of this plasma is very low 
so that the mean free path for collisions is much longer 
than the dimensions of the obstacle itself, the earth's 
magnetosphere. Therefore a classical collisional shock 
wave cannot exist. Satellite observations showed the 
presence of a well-defined shock transition of the order 
ei too Km. 

In the theory of shock waves, it is important to determine 
the properties of the medium after passage of the shock 
in terms of the unperturbed quantities before the shock 
Such as densities, pressure and temperature. Assuming that 
Ene Shock profile is (quasi-)stationary in a certain 
reference frame, a number of relations between upstream 
and downstream quantities are given by jump conditions, 
which derive directly from the conservation of mass, momen- 
tum and energy. In usual gas dynamics and MHD-theory, 
these jump conditions are sufficient to determine the basic 
properties of the state behind the shock although a theory 
Of the shock structure and the dissipation process leading 
Bementropy production requires a kinetic description. 

In case of a collision-free shock wave the plasma is, 
in general, not in local thermal equilibrium behind the 
shock transition. Electron and ion pressures may be quite 


different from each other, depending on the dissipation 
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mechanism. Here, the jump conditions do not completely 
describe the state behind the shock. A dissipation of the 
Shock structure itself being needed. This would, e.g. 
provide information on the ion or electron thermalization 
length. 

Plasma momentum transfer calculated to be collision 
free was observed only when a magnetic field was present 
by Friedman and Dawson [Ref. 11] during plasma gun studies. 

Dean et. al. [Ref. 8] performed an experiment which was 
Similar to this present one. A Nd-doped glass laser 
(8J, 30 nsec) was used. The target was 0.25 mm Lucite 
(C,H,0.) fiber suspended in an ambient nitrogen gas. They 
found collisionless interactions even in the absence of 
an applied field and explained this apparent discrepancy 
by the fact that spontaneously generated fields always 
existed in their experiment. They meaSured the thickness 
of the interaction front as approximately 1mm and esti- 
mated the binary-collision momentum transfer mean free 
path, from a screened coulomb potential of the order of a 
centimeter. 

Shock thicknesses in collisional plasmas are usually 
taken as one-to-several collisional mean free paths. The 
difficulty of establishing a given experimental situation 
as collision free has been emphasized by Wright [Ref. 32] 
who suggested that the results of [Ref. 8] can be explained 


by collisional coupling alone. Dean et. al. [Ref. 9], 
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after Wright's suggestion reperformed their experiments. 
They changed the target to pure carbon filaments to find 
out whether the multicomponent nature of the former target 
(C.H,0,) was playing a role. As a result they found the H 
and O in the target was inconsequential. They also changed 
background gas and used Helium instead of Nitrogen. The 
basic phenomena were qualitatively unchanged. They sug- 
gested that the momentum coupling region which was observed, 
was not a shock. It was merely a region of momentum 
coupling between interstreaming plasmas. The mean free 
path 1s the e-folding distance for momentum coupling, so that 
several mean free paths are required for ahighdegree of coupling. 
The expanding olasma heats and transfers momentum 
Bemeene photoionized background plasma. These inter- 
actions can be either collisional or collision-free depending 
On densities and velocities. During an experiment per- 
formed by Stamper et. al. [Ref. 27], the background densi- 
ties have typically been such that electron processes are 
collisional while the ion processes are collision-free. They 
Observed a sharply defined momentum coupling interaction 
region between the laser plasma and background plasma, but 
they have not observed the existence of a shock. 
Koopman and Tidman [Ref. 17] observed a thin collision- 
less electrostatic shock front propagating radially at the 
leading edge of the plasma ball. In their investigation, 


=>granules of LiD have been irradiated by 15-nsec, 3-J pulses 
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from a Ruby laser. The target consisted of a cluster of 
5-to-25 um diameter solid particles which were injected 
into the 400-um diameter laser focal region by an electro- 
mechanical apparatus. No magnetic fields were present. 
The phenomenon is analogous to a blast wave; the expanding 
plasma acts like a piston and sweeps up the background 
plasma. This swept-up plasma gives rise to a spherical 
lon-wave shock. 

The distinction between collision-free and collisional 
mechanisms in counterstreaming plasma studies, and the 
appropriate techniques for calculating collisional cross 
sections for momentum transfer, have been the subject of 
considerable controversy. The question of the role of 
binary collisions is always a very complicated one in plasma 
experiments. There is still a need for more sophisticated 


treatments of collisional effects in plasmas. 
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ie EXPERIMENTAL METHODS 


A. LASER 

The laser used in this investigation was a Korad-K-1500 
O-switched neodynium-doped glass laser. The laser consisted 
of an oscillator and a single amplifier stage. The laser 
was Q-switched with a Pockels cell. The laser system was 
Operated at a power level of 230 MW, which was obtained 
with a pulse width (full width at half-maximum) of 22 nano- 
seconds. The output of the system varied within about 10 
percent of 230 MW and no data were accepted if obtained at 
power levels outside this range. Other details of the laser 
are given by Davis [Ref. 6]. The laser focal spot size is 
approximately 0.15 cm in diameter after focuSing with a 


Eeeem focal length lens. 


SeeeeuaAskR MONITORING TECHNIQUES 

The laser beam was monitored with a Korad-K-Dl photo- 
@tede. This photodiode provides a signal which is propor- 
tional to the laser power and also by integrating, gives 
another signal which iS proportional to the laser pulse 
energy. 

About four percent of the beam was reflected by a beam 
pemeeter to a MgO diffusing block for viewing by the photo- 
diode. The scattered light then passed through a 0.1 per- 


cent neutral density filter before entering the photodiode 
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as shown in Figure 2. The power Signal was used to trigger 
a Tektronix 555 Dual-beam oscilloscope and was also dis- 
played on that oscilloscope. 

The energy Signal was displayed on a Tektronix 564B 
storage oscilloscope. 

The photodiode energy Signal was calibrated by L.L. 
McKee [Ref. 18] by uSing a Westinghouse RN-l laser radiometer. 
The radiometer provided an absolute measure of the energy 
incident at the target which allowed calibration of the 
diode's integrated signal. 

The energy output of the laser was monitored on every 
SmOe.) FOr energy outputs within 10 percent of 5.0 Joule, 
the pulse width was consistently 22 + 1 nsec. A two minutes 


cooling period between shots resulted in good reproducibility. 


eee EMPORAL REFERENCE FRAME 

The temporal relation between the diagnostic probe 
response and the leading edge of the laser pulse was deter- 
Mined by taking into account the various cable delays and 
laser beam optical path lengths. A long (10 meter) delay 
line was used from the diagnostic probe to the oscillo- 
Scope in order to separate the laser power signal and early 
diagnostic responses which would otherwise be superimposed 
On each other. 


10 cm/sec 


By asSuming that the Signal travels at 1.98 x10 
in the RG-174/U coaxial cable, the delay was computed as 


34 nsec. 
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D. PLASMA CHAMBER 

The target was located in a chamber which was specially 
constructed by McKee to be used for the investigation of 
laser-plasma. The inside volume of the chamber was cylin- 
drical with a diameter of 10 inches and a height of 6 
inches (Figure 3). 

The beam entered the chamber after passing through a 
30 cm focal length lens and hit the target at an angle of 
30° so that one can probe the plasma along the target nor- 
mal which appears to be the axis of greatest symmetry of 
the expanding plasma. 

The target used was an aluminum flat disc, 1/4 inches 
thick, and about 2 inches in diameter. It was rotated after 
a-proximately 10 shots to prevent severe cratering of the 
target in the laser impact area. 

The chamber was pumped to a vacuum of about pian” Owe 
of air and then background gas was allowed in by a needle 
valve to fill the chamber up to the deSired pressure. The 
Meecuum Of about eens Torr was obtained by using an oil 
diffusion pump. Vacuum pressures lower than 0.1m Torr 
were measured with an lonization gauge. A thermocouple 
gauge was used to measure pressures between 0.1m Torr and 
1000 m Torr, while, above 1 Torr pressures were monitored 


by a manometer filled with diffusion pump oil. 
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BE. MAGNETIC PROBES 
ie) ProbesGalibration 

The magnetic probe used in this investigation con- 
sisted of five turns of #40 Formvar copper wire. The probe 
coil was about 1mm in diameter. The probe circuit response 
was linear, with an effective area (nA) of the inductive 
ieom Of 517 LO o ma? for frequencies below about 8 MHz. 
The effective area (nA) of the probe was measured by using 
Helmholtz coils. The instruments employed for this measure- 
ment were a General Radio 1330A continuous wave oscillator, 
Hewlett-Packard 411A r.f. VTVM, Helmholtz coils and a 
Tektronix 545B oscilloscope. A circuit diagram for probe 
calibration is shown in Figure 4. 

The details of construction of Helmholtz coils and 
the use of these coils for calibration of a probe is des- 


cribed by McLaughlin [Ref. 19]. The effective probe area, 


according to McLaughlin is given by the equation 


where K is a parameter which differs only from coil to 


coil (for the coil which was used in this investigation, 


me 0.025 m-), Eis the voltage across the probe connection 
cable and Vr is the voltage applied to the Helmholtz coil. 
2. Probe Signal Reliability and Sensitivity 


To insure that the magnetic probe signals (dB/3t) 


were meaningful, several steps were taken. The probe 


47 





Signals were monitored under experimental conditions for 
which no target was in place, to verify that a null signal 
was obtained. Also, the probe was rotated by 180° about 
its axis and the probe signals were checked to make sure 
that they reverted in polarity. This demonstrated that 
the signals were not electrostatic in nature. Magnetic 
field signals were recorded on a Tektronix 555 dual-beam 
oscilloscope. 
3. Probe Signal Integration 
Pecoraung —EO,eleceromagnectre theory, the output of 


a@coil is 


eke 


where nA 1s the effective coil area, B is the magnetic 
field, and e(t) is the voltage generated across the coil 


which was recorded on a Tektronix 555 dual-beam oscillo- 


scope. Integration of equation (26) gives: 
B t 
= a 
B = f dB = =F ieee) at (27) 
0 0 
1: 


The integration of the probe signal ( { e(t)dt) was 
0 


performed by meaSuring the area under the e(t) curve. 
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foes lect RiC DOUBLE PROBES 

Double probe thoery can be found in several references 
Weent. 3) and [Ref. 4]. 

At sufficiently high bias voltage, the probe draws 


mee Lon Saturation current I 


I ~ A, Ns Ze Vip (28) 


where A, is the probe area perpendicular to the direction 


of plasma flow, N.Ze is the average ionic charge per unit 


volume in the laser plasma and V 1s the magnitude of the 


Lp 


plasma flow velocity. 


In equation (28) the assumptions are: 


(1) The areas of probe tips are equal. 
¢ 
(2) The plasma parameters are uniform across the probe 
eps. 
(3) Plasma-probe interactions such as secondary elec- 


tron emission can be neglected. 

The secondary electron emission coefficient (¢€) due to 
10on bombardment (Y.), defined as the average number of 
electrons ejected per incident ion, has been measured by 
Koopman [Ref. 16] for similar experimental conditions and 
meena tO fall in the range of 0.0 - 0.5. Therefore secon- 
dary electron emission will contribute to the total current 
drawn between the probe tips. After taking secondary 


electron emission into account equation (28) becomes 
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eae ie Ue) Ay n; Ze V 2) 


Lp 


The probe was biased at 20 volts to collect a saturated 


ion current. A large secondary emission coefficient of 


¢ = 0.5 was used in obtaining a quantitative idea of plasma 
density. 
The probe circuit 1S Shown in Figure (5). The probe 


circuitry was enclosed by an aluminum housSing attached to 
the probe and the entire system was isolated from ground. 
The diameter of the cylindrical tips was about 1 mil, they 
had a length of 1 mm and were about 1 mm apart. 

The double probe signals were recorded on a Tektronix 
555 dual-beam oscilloscope. To obtain an estimate of the 


plasma density, equation (t) 1s solved for Ns and evaluated 


with Z = io, © =40.5, Ay = DBA se MS oa, I= = Pee ee 
a 22 V (volts) 
My = 1209 x10 x v Pen eES) (30) 


where V is the voltage across a ] ohm~-resistor (see Figure 5). 


fee DaArtA COLLECTION AND ERROR ESTIMATES 

The magnetic probe was inserted into the vacuum chamber 
through a small pipe which was soldered onto a brass base 
plate. 

The base plate could be moved about on an "O" ring seal 
located on the top of the chamber without breaking the 


vacuum seal. The base plate was scribed with reference lines. 
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These lines allowed positioning of the plate over a one 
millimeter grid attached to the top of the chamber. The 
position of the probe relative to the target was measured 
internally at the beginning of each data session to insure 
that the probe had not moved in its holder. This method 
allowed positioning of the probe to within + = mm of the 
stated position in relation to the target. The electric 
probe was inserted into the chamber from a side window. 

The relative measure of plasma density from point to 
point 1S considered good due to the high reproducibility 
of the laser pulses. It is eStimated that these relative 
values are accurate to within +10 percent. 

The laser gave very reproducible pulses; the output 
energy was allowed to vary as much as + 10 percent. 

The area of probe tips is considered accurate to within 


10 percent. The plasma expansion velocity is considered 


accurate to within 10 percent. 
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V. EXPERIMENTAL RESULTS 


During this investigation, hydrogen gas was used as 
the background gas. An electric probe was introduced from 
the port No. 6 (Figure 3) and the magnetic probe from the 
top window. 

To decelerate the plasma front rapidly, a glass plate 
1s placed into the vacuum chamber (Figure 7). The main 
influence of the glass plate on spontaneously produced 
magnetic field is displayed in Figure 28. This figure 
indicates that the maximum value of reverse field when the 
glass plate is in place (at z = 1.15 cm), measured with 
the probe at r = 0.4 cm, 9 = 0° and z = 1.0 cm, is greater 
than the maximum value of reverse field when the glass 
mere 1S not in place (at Zz = 7.5 cm). 

Figures 8 and 9 display the relation of sB/3st to the 
plasma density for a vacuum pressure of aon Torr, as 
a function of time. When the glass plate is in place, the 
Magnitude of dB/dt is greater than the magnitude of dB/ot 
when the glass plate is not in place. 

Figure 16 demonstrates that the slope of the plasma 
density curve when the glass plate is at z = 1.15 cm is not 
very different from the slope of the plasma density curve 
when the glass plate 1s at z = 7.5 cm 

Figure 17 shows that when the background gas pressure 


fw. 5 Torr (H5), the existence of the glass plate at 
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z = 1.15 cm causes a noticeable steepening and a sudden 
increase in plasma density. 

Figures 10 and 11 show the relation of 3B/dt to the 
plasma density at a background pressure of 0.5 Torr (H.) 
as a function of time for two locations of the glass plate. 

When the background gas pressure is 5.0 Torr (H,), the 
Sudden increase in plasma density when the glass plate is 
in place (Figure 18) is much more pronounced (see also 
Beeguires 13, 15 and 28). For 5.0 Torr (H,) background 
pressure, the difference between the maximum value of 
reverse field when the glass plate is at z = 1.15 cm and the 
maximum value of reverse field when the glass plate is 
at z = 7.5 cm is greater than the difference between the 
Same quantities for a background pressure of 0.5 Torr (H,). 

Figures 12 and 14 display the relation of dB/dt to the 
plasma density for 5.0 Torr background gas pressure at the 
eommeer = 0.4 cm, 9 = 0° and z = 1.0 cm for two different 
locational of the glass plate (note that each curve in this 
Figure was multiplied by a different factor). 

Figures 19, 21 and 23 demonstrate the comparison of two 
plasma density curves each for different background gas 
Smeessure, at the point r = 0.4 cm, 9 = 0° and z = 1.0 cm 


with the glass plate at z = 1.15 cm. The slope of each of 


these two curves is almost the same. The comparisons of 
the same curves when the glass plate is at z = 7.5 cm are 
displayed in Figures 20, 22 and 24. From these figures, one 
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can notice that the slopes of the two density curves each 
one for different background gas pressures are different 
from each other. 

Figures 25, 26 and 27 illustrate the comparison of 
the magnetic probe Signals (dB/dt) when the glass plate 
1s located at z = 1.15 cm with the magnetic probe Signals 
(oB/ot) when the glass plate is located at z2 = 7.5 cm for 
2x10 ° Torr, 0.5 Torr (H,) and 5.0 Torr (H,) background 
pressures. From these figures, one can see that the reverse 
field maximum is the greatest when the background gas 
pressure is 5.0 Torr (H.) and glass plate is at z = 1.15 cm. 

Figures 29, 31 and 33 illustrate the comparison of the 
two magnetic probe Signals (dB/dt), each one for a differ- 
ent background pressure, as a function of time at the point 
fees cm, € = OY and z= 1.0 cm and the glass See is 
in place. An investigation of these figures indicates that 
the higher the background gas pressure, the greater the 
area under the negative portion of the dB/dt curve, while 
the area under the positive portion of 3B/dt curve becomes 
smaller. This means that as the background gas pressure 
1S increased, the reverse field maximum becomes greater when 
the glass plate 1s in place. Figures 30, 32 and 34 illus- 
trate the same comparison when the glass plate 1s at 
ge /.5 CM. 

For a vacuum pressure of 2x10 ° Torr, the magnitude 


of the magnetic probe signal (dB/dt) when the glass plate 


1S in place is about 4 times greater than the magnitude of 
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the magnetic probe Signal (dB/dt) when the glass plate is 
not in place (see Figure 25). (Note that each curve in 
this figure was multiplied by a different factor.) 

When the glass plate is in place, the magnitude of the 
Magnetic probe Signal (dB/dt) at a vacuum pressure of on 
Torr is greater than the magnitude of the magnetic probe 
Smonals (dB/dot) at 0.5 Torr (H.) background pressure and 
fies .0 Torr (H.) background pressure (see Figures 29 and 
31). But when the glass plate 1S not in place, the magnitude 
of the magnetic probe Signal (dB/3dt) at Dee ® Memmi aviee Ulm 
pressure is smaller than the magnitude of the magnetic 
probe Signals (dB/dt) at 0.5 Torr (H.) background pressure 
mceat 5.0 Torr (H.) background pressure (see Figures 30 
and 32). From Figure 34, it can be seen that when the glass 
plate is not in place, the magnitude of the magnetic probe 
eeemal (oB/dt) at 0.5 Torr (H,) background pressure is 
greater than the magnitude of the magnetic probe signal 
mMeeyvot) at 5.0 Torr (H.) background pressure. (Note that 
each curve in this figure was multiplied by a different 
factor.) But Figure 33 illustrates that when the glass 
plate is in place the magnitude of the magnetic probe signal 
(9B/dt) at a background pressure of 0.5 Torr (H), between 
0 and 100 nsec is greater but between 100 and 450 nsec it 
is smaller than the magnitude of the magnetic probe signal 
Mee7ot) at a background pressure of 5.0 Torr (H.,) . (Note 
that each curve in this figure was multiplied by a different 


jmoaeccor. ) 
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Figure 35 illustrates that when the glass plate is 
in place, for a background pressure of 5.0 Torr (H5), the 
Slope of the plasma density curve at the point r = 0.4 cn, 
89 = 0° and z = 1.0 cm iS greater than the slope of the 
plasma density curve at the point r = 0.8 cm, 6 = 0° and 
z= 1.0cm. But when the glass plate is not in place 
Figure 36) the slopes of two plasma density curves are the 
same. 

During this investigation, for a background pressure of 
ero) LOrr (H5), the maximum value of the magnitude of the 
reverse magnetic field was measured as 150 gauss when the 
glass plate is in place at the point r = 0.4 cm, 8 = 0° 
and z = 1.0 cm and at t = 440 nsec, while at same background 
pressure, point and time but when the glass plate 1s not 
in place, the maximum value of the magnitude of the reverse 


magnetic field was measured as 45 gauss only. 
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Vis COMPARISON OF THEORY AND EXPERIMENTAL RESULTS 


SoemenieLonead in Secelom FII-C, the production of initial 
field requires a density gradient in the -zg direction and 
a temperature gradient in the-rdirection (see Figure 6). 
The cross product of these two vectors generates a magnetic 
Pe—meein the © direction 1.e., In the clockwise direction 
when looking in the direction of the laser beam. 

If one of the gradients 1s reversed tothe opposite 
direction, a magnetic field, in the counterclockwise 
direction is produced. The radial temperature gradient is 
always in the -r direction and its direction remains unchanged, 
either the glass plate is in place or not. Therefore the 


only possible gradient which can be reversed by placing 


the glass plate at z = 1.15 cm, is the axial density gradi- 
emee ©€he equation of continuity is 
ony - (nv) = 0 10 
ot 


if cylindrical symmetry is assumed; 
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Also for the laser-produced plasma expanding in the z- 


emeection 





--. a on iy where yemsmthe expansion velocity of 


thus: is 
the plasma front. 

PisieemMmcenomsriateamrene L£elation Of Be to the 
plasma density as a function of time at the point r = 0.4 cm, 


8 = 0° and z = 1.0 cm for a background gas pressure of 


eo TOrr (H.) when the glass plate is at z = 1.15 cm. 
dn 


Between t =~ 100 nsec and t ~ 180 nsec. es 1S positive, 
which means that “= 1s negative or the axial density 


gradient is in the-z direction. Figure 28 illustrates 
that at t = 150 nsec the reverse field appears when the 
glass plate is in place. Since the axial density gradient 
at this time is in the-z direction, the production of the 
reverse field can not be correlated to any change of Sign 
in the direction of the axial density gradient, since 
there is no such change when the reverse field is "turned 
Om. 

There always exists a radial density gradient in the -r 
direction, no matter where the glass plate is (see Figures 
35 and 36). Therefore it requires an axial temperature 
gradient in the -z direction. Such an axial temperature 
gradient can be generated by the shock heating of the 
plasma. 

Figure 28 demonstrates that the magnitude of the 
initial magnetic field, when the glass plate is in place 
is smaller than the magnitude of the initial magnetic 


field when the glass plate is not in place. This decrease 
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in magnitude of the initial magnetic field is caused by 

the decrease in axial density gradient due to pile-up of 
laser plasma in front of the glass plate. This piled 

up plasma then diffuses radially and also causes a decrease 
in radial density gradient. This decrease in radial den- 
Sity gradient can be seen if one compares Figure 35 and 
Pergure 36. 

Figure 18 displays the steepening in the plasma den- 
Sity between 150 and 200 nsec, when the expanding plasma 
impinges on a glass plate. This steepening can be related 
to the formation of a shock wave. Such a shock wave causes 
an increase in plasma temperature (see Appendix A) and this 
increase results in an axial temperature gradient in the 
~z direction, since the shock wave travels in the +z direction. 


As long as the ratio, 


OB 
sl aisles 1 
eee ee | 
ot 


1s greater than one, the reverse field only reduces the 
magnitude of the initial magnetic field. But when this 
ratio decreases to a value which is less than one, the 
reverse field becomes the dominant field. 

From Figure 18, the ratio of two densities, one at 


feet a2.5/ nsec the other at t = 180.83 nsec, is: 


ee SS and At = 21.46 nsec 
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in order to obtain this density ratio, the shock wave 
must have a Mach number of the order of 1.37. According 
to the classification of Jaffrin and Probstein [Ref. 15] 
this shock is a "moderate" shock whose thickness is of 
the order of %/e, where 2 is the downstream (before shock) 
—€) 1/2, 


ion-ion collision mean free path, ande = a The 


electron temperature increases in this ever the ion 
temperature increases in the Second inner layer (see Section 
mero E) whose approximate thickness is the downstream ion- 
10n collision mean free path 2. 
For this Mach number M = 1.37, a temperature ratio, 
ve 
ve 
is obtained from the Rankine-Hugoniot jump conditions 
(see Appendix A, equation 3). 

Since the shock thickness can be defined from the 
electric probe data as 56 = Vt, where V is the velocity of 
the expanding front and t is the rise time of the probe 
Signal. 

If no interactions were occurring between the inter- 
streaming plasmas, then the plasma expansion velocity 
would be independent of background gas pressure. The 
Observed dependence of expansion velocity can be seen in 
Figures 20, 22, 24, between 100 and 200 nsec. From these 
figures one can conclude that as background gas pressure 


increases, the expansion velocity decreases. 
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To have an idea about the structure of the shock 
one must compare the shock thickness with the collison 
mean free path length. The ion-ion collision mean free 


path can be calculated from 


Ig 
2 
(en) = want xT. (eV), 
n, (cm <n 
where T. is the ion temperature. In order to calculate 2, 


One needs to know ns and T; - 

Boom kT, = 5m.V*, where V is the expansion velocity 
and m 1s the ion mass, T; can be estimated. If one 
assumes that the Shock wave propagates in the background 
gas, or background gas flows into the shock relative to 
the frame of the shock, for hydrogen background gas, 

a 
m. =m = 1.673 x ou gG and V = 10/ cm/sec. For these 
values of m5 and V, T; = 52 eV is obtained. Assuming a 
fully dissociated and fully lonized background gas at 
5.0 Torr the downstream ion density 1s 
16 =e ey ~3- 


ae eo en LO cm oS a 0 cm 


For these values of ne and Tis {*1s Of Ehe order of 
8.2 x 107°? cm. 

The shock thickness can be obtained from Figure 18, 
in which t = At = 21.46 x ig sec and for V ~ Tron cm/sec 


Pees TV = 21.46 x 107° — 0) 2 Mek ihe Since the shock thickness 6 
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is greater than the ion-ion collision mean free path, 

the Povametesvann of a collisional shock wave must be expected. 
In summary it can be said that whenever the background 

gas pressure is high enough, the expanding laser plasma 

acts like a piston and sweeps up the photoionized back- 

ground plasma. In front of this piston a shock is formed 

to drive an axial temperature gradient which is large 


enough to produce a reverse magnetic field. 





VII. SUMMARY AND CONCLUSIONS 


The observed dependence of the plasma front expansion 
velocity on background gas pressure indicates the exis- 
tence of an interaction between the laser plasma and the 
photoionized background plasma. As a result of this inter- 
action, a collisional shock wave, driven by the laser 
plasma piston is formed. This shock wave heats and com- 
presses the plasma. 

Initially the spontaneously produced magnetic field is 
in the §-direction (see Figure 1) and the thermal source 
term for these magnetic fields consists of the cross 
product of the radial density gradient and the axial den- 
Sity gradient. When the shock wave is formed an axial 
temperature gradient is produced. This axial temperature 
gradient and the radial density gradient which always 
exists in the Silom Shaped laser plasma are the main cause 
of the reversal of these initial magnetic fields. 

The comparison of the shock thickness and the colli- 
sion mean free path determines the type of the shock. If 
the shock thickness is greater than the collision mean 
free path, the shock is collisional in nature. 

Determining the mean free path length in a plasma is 
always a complicated undertaking and the results should be 
considered as crude estimates. More sophisticated treat- 
ments are necessary in order to determine the collisional 


effects in plasmas. 
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The reversal of the spontaneous magnetic field indi- 
cates the existence of a shock wave propagating through 
the plasma. Such generation of magnetic fields by shock 


waves can be a very important mechanism in astrophysics. 
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VIII. RECOMMENDATION FOR FURTHER STUDY 


The electron temperature has been evaluated indirectly 
from the shock parameters. An experimental measurement 
of the electron temperature gradient is required for a 
further study of field reversal in shock waves. Plasma 
temperature, before and after the plasma front encounters 
the glass plate can be meaSured by means of x-ray diagnostic 
techniques or by other means. Also, more accurate data 
about the r dependence of the plasma density are needed 
to provide a quantitative comparison between experiment 
and theory. The data about the r dependence of the plasma 
density can be obtained by use of the electric probe 


technique. 
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APPENDIX A 


Rankine-Hugoniot jump conditions can be determined by 
uSing the equation of continuity, conservation of momentum 
density and conservation of energy density in the plasma 


for steady state conditions [Ref. 29]. 
S 
The Mach number is defined for = = yY = 2 as; 
V 


Ws 
M = - 


Pail 
3M 


) 
] 


where Vi 1s the downstream (before shock) fluid velcoity, 


Te 1s the downstream electron temperature, Tod 


downstream ion temperature and M is the ion mass. 


1s the 


foe UMP COndic10ns are: 


n V 
a 2 ae 
cae 7 oo (2) 
N>5 Vy 4 us 
Bag Easy 8 (Nema (Ee : 
T +7 a a 
il “el L6M 


where the subscript "2" implies upstream (after shock) 


value. 
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Figure 2. Block diagram of experimental layout 
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meogure 3. 





Top view of vacuum chamber. Port #1 is the laser 
beam entry port, #2 1s the reflected laser beam 
observation port, #3 is the transmitted laser 
beam observation port, #4 is the optical observa- 
ttOneport, #2 2s the rglass plate entry port, 

#6 is the electric probe entry port, #7 is the 
probe observation port. 
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Figure 6. Geometry of a laser produced plasma and 
the self generated magnetic field 
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